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Ab  initio  quantum  chemical  calculations  have  been  performed  to  study  the  adsorption  of  trichlorophosphate, 
dimethyl  methylphosphonate,  and  Sarin  via  hydrogen  bonding  to  Si— OH  groups  on  the  amorphous  Si02 
(rz-SiCL)  surface.  Two  Si03  models  are  used:  a  small  SLChHx  “cagelike”  cluster  and  a  larger  SLiOmPLx 
structure  designed  to  approximate  the  local  environment  in  a-SiCK  In  the  latter  case,  regions  of  different 
local  OH  density  are  considered  as  adsorption  sites.  Adsorption  energies,  bonding  geometries,  and  adsorbate 
vibrational  modes  are  obtained,  and  anharmonicity  is  explicitly  included  in  the  treatment  of  the  SiO— H 
stretching  mode.  The  computed  results  for  the  adsorption-induced  shift  in  frequency  of  the  SiO— H  stretch 
and  of  the  molecular  P=0  stretch  are  compared  with  the  available  experimental  data.  For  all  three  species, 
the  most  stable  adsorption  geometry  involves  hydrogen  bonding  between  two  Si— OH  groups  and  the  O  atom 
of  the  P=0  group. 


1.  Introduction 

The  interaction  of  chemical  warfare  agents  (CWAs)  with 
materials  is  an  important  issue  in  protection  and  remediation. 
Of  specific  concern  is  a  basic  understanding  of  the  long-term 
stability  of  such  species  (termed  “agent  fate”)  under  environ¬ 
mental  conditions  and  of  the  nature  of  the  chemisorption  bond 
at  the  material  surface.  These  subjects  are  difficult  to  address 
experimentally  because  of  the  extreme  hazards  in  working  with 
real  CWAs.  Relatively  benign  simulants  are  typically  used 
instead,  but  it  is  often  uncertain  the  extent  to  which  such  species 
mimic  the  adsorption  properties  of  real  CWAs.  A  direct 
evaluation  of  a  particular  simulant  involves,  again,  experimenta¬ 
tion  with  the  corresponding  real  agent. 

An  obvious  approach  to  this  problem  is  to  apply  quantum 
chemical  (QC)  methods  to  model  the  systems  of  interest.  This 
has  been  done  in  a  number  of  cases1-11  involving  mainly  the 
interaction  of  CWAs  with  ionic  oxides  or  with  molecular 
species.  For  most  metal  oxides,  adsorption  generally  involves 
Lewis  acid— base  chemistry  for  OH-free  surfaces  and  Brpnsted 
acid— base  chemistry  for  hydroxylated  surfaces.  On  the  other 
hand,  for  Si02  and  related  materials  hydrogen  bonding  (H 
bonding)  to  surface  silanol  (Si— OH)  groups  is  involved,  and 
accurate  theoretical  treatment  of  these  relatively  weak  surface 
H  bonds  can  be  difficult.12  Recently  the  adsorption,  via  H 
bonding,  of  CWAs  on  aluminosilicate  clay  surfaces  has  been 
modeled.3,9  The  structures  formed  are  complex  and  depend  on 
the  extent  of  hydroxylation  of  the  surface. 

The  present  work  applies  ab  initio  QC  methods  to  the 
interaction  of  three  related  phosphate  and  phosphonate  molecules 
on  amorphous  silica  (a-SiCh).  These  are  (Figure  1)  trichloro¬ 
phosphate  [Cl3P=0,  TCP],  dimethyl  methylphosphonate  [(CH30)2- 
(CH3)P=0,  DMMP]  and  Sarin  [(iPrO)(CH3)(F)P=0,  where  iPr 
=  (CH3)2CH-].  TCP  is  used  here  as  a  relatively  simple  species 
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Figure  1.  Schematic  diagram  of  the  molecular  structures  of  TCP, 
DMMP,  and  Sarin. 


for  testing  the  computational  procedures  and  for  probing  the 
reactivity  of  various  parts  of  the  model  clusters.  DMMP  is 
important  as  a  safe  simulant  for  the  nerve  agent  Sarin.  The  main 
focus  of  this  work  is  in  comparing  theoretical  results  for  the 
adsorption  of  TCP  and  DMMP  on  a-SiCL  with  experimental 
data,  especially  the  infrared  (IR)  vibrational  spectrum  and  in 
comparing  results  for  the  adsorption  of  DMMP  versus  Sarin. 
The  goals  are  an  evaluation  of  the  reliability  of  QC  methods  in 
describing  the  adsorption  and  also  a  “side-by-side”  comparison 
of  the  behavior  of  the  simulant  and  the  real  agent. 

2.  Computational  Methods 

All  calculations  were  done  using  the  Gaussian03  suite  of 
programs13  with  standard  basis  sets  built  into  the  package.  Most 
of  the  work  was  done  with  6-31+G(2d,2p)  basis  sets,  but  some 
results  were  obtained  for  comparison  at  the  6-31 1++G(d,  p)  and 
6-31+G(d,  p)  levels.  It  is  clear  from  previous  work14,15  that 
heavy-atom  diffuse  functions,  together  with  polarization  func¬ 
tions  both  on  heavy  atoms  and  on  H,  are  necessary  for  an 
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adequate  description  of  H  bonding  in  the  present  case.  Hence, 
6-31+G(d,  p)  is  considered  to  be  the  smallest  basis  set  that  can 
be  used  reliably.  It  has  also  been  shown15  in  studies  of  NH3 
adsorption  on  SiCT  that  6-31+G(2d,2p)  basis  sets  yield  a 
somewhat  more  accurate  adsorption  energy  (A£ads,  see  below) 
than  do  6-31+G(d,  p)  but  that  further  enlargement,  up  to 
6-31 1++G(2d,2p),  has  no  significant  additional  effect  on  AE^. 
All  calculations  were  done  using  density  functional  theory  (DFT) 
with  the  B3LYP  functional.  In  the  past,  questions  have 
sometimes  been  raised  concerning  the  applicability  of  DFT  to 
H  bonding.  However,  several  previous  studies  (e.g.,  refs  14— 
20)  have  demonstrated  the  validity  of  DFT  with  the  B3LYP 
functional  in  describing  the  H  bonding  of  adsorbates  (including 
phosphonates1419)  on  Si02  surfaces. 

Much  of  the  present  work  involves  a  comparison  of  observed 
and  calculated  IR  vibrational  spectra  as  a  means  of  evaluating 
the  computational  methods  and  results.  Particular  attention  is 
given  to  the  shift  to  lower  frequency,  Ai'(0— H),  of  the  SiO— H 
stretch  that  results  from  adsorption.  Previous  work15’16-21  has 
shown  that  accurate  results  for  Av(0— H)  require  that  anhar- 
monicity  be  treated  explicitly.  Computation  of  anharmonic  terms 
for  all  3N-6  modes  of  an  N-atom  model  (adsorbate  plus  SKH 
cluster)  is  not  feasible  for  systems  of  the  size  studied  here; 
therefore,  the  approach  described  previously22  is  followed.  The 
v(0~ H)  mode  is  well  separated  in  energy  from  (and  therefore 
decoupled  from)  all  other  modes.  Beginning  with  the  relaxed 
structure,  a  series  of  single-point  calculations  is  done  in  which 
the  H  atom  is  displaced  incrementally  in  a  direction  parallel  to 
the  O— H  bond.  In  total,  28  displacements  in  the  range  of 
—0.55  <  (5o-h  —  +0.85  A,  relative  to  the  equilibrium  O— H 
distance  of  about  0.96  A,  were  used. 

The  one-dimensional  (ID)  energy  surface  thus  obtained  was 
least-squares  fitted  with  a  polynomial  of  order  12  <  k  <  15  in 
(5oh,  and  the  resulting  potential  function  was  then  used  to  solve 
the  ID  Schrodinger  equation  numerically.23  A  high-/:  polynomial 
was  necessary  to  optimize  the  fit  over  a  wide  range  of  (5o-h-  A 
similar  procedure  was  also  used  to  estimate  the  effects  of 
anharmonicity  in  the  P=0  stretching  mode  of  TCP.  The  mass 
used  for  the  oscillator  is  the  reduced  mass,  /uxy  =  (inxmy'lKmx 
+  my).  From  the  lowest  three  eigenvalues  (Eo,  Ei,  and  E2)  the 
quantities  of  interest  are  obtained.  These  are  the  anharmonic 
fundamental  and  first-overtone  frequencies  voi  =  (Ei  —  Eo)  and 
V02  =  (E2  —  Eo)  and  the  anharmonicity  parameter  ryte  =  (2v0i 
—  v02)/2  where  ve  is  the  harmonic  frequency. 

The  approach  described  above  involves  an  approximation  that 
should  be  valid  in  the  case  of  v(0~ H).  In  obtaining  the  ID 
potential  energy  surface,  the  bond  stretching  is  modeled  entirely 
in  terms  of  H-atom  displacement;  whereas,  the  true  normal  mode 
is  more  accurately  described  in  terms  of  the  individual  atomic 
displacements  as  (5o-h  =  0.94(5h  —  0.06(5o-  Different  effects 
on  the  total  energy  are  expected  for  equivalent  displacements 
of  H  and  O  ((5h  and  <5o),  but  any  error  in  the  ID  potential 
function  resulting  from  the  approximation  (5o-h  ^  (5h  has  been 
neglected. 

The  frequencies  of  other  modes  are  obtained  in  the  harmonic 
approximation  and  scaled  using  factors  given  in  previous  work. 
The  scaling  accounts  for  anharmonicity  in  an  approximate 
manner,  based  on  an  average  correction  for  a  wide  range  of 
molecules  and  normal  modes.  For  B3LYP/6-31+G(d,  p)  the 
factor24  is  0.9632,  and  for  B3LYP/6-311++G(d,  p)  the  factor25 
is  0.9679.  To  our  knowledge,  no  scaling  factor  has  been  reported 
specifically  for  B3LYP/6-31+G(2d,2p);  hence  the  factor  of 
0.9632  for  B3LYP/6-31+G(d,  p)  will  be  used.  The  uncertainty 
in  these  factors  is  in  the  range  of  ±0.02—0.03  (ref  24). 


(a)  (b) 


Figure  2.  Schematic  diagram  of  the  model  a- SiO?  cluster.  A  model 
reactive  region  with  a  relatively  high  OH  density  is  shown  in  (a).  The 
same  reactive  region  is  shown  in  (b)  as  part  of  the  total  cluster.  Si  (O) 
atoms  are  shown  in  blue  (red)  except  for  the  Si  atoms  at  the  reactive 
region  which  are  in  green.  H  atoms  are  gray.  The  structure  shown  in 

(a)  constitutes  the  ONIOM  model  system,  and  the  structure  shown  in 

(b)  is  the  ONIOM  real  system.  Further  details  are  given  in  the  text. 

Other  results  of  interest  are  A£ads  and  related  quantities, 
defined  by 


A/-’ac|s  =  TYmol  +  cluster)  —  TYmol )  — 


A^ads  A^ads  ' 


A  Ea 


A^ads  A£ads+  A-EzPE 


^(cluster) 


where  the  first  term  in  A£ads  is  the  total  energy  of  the  relaxed 
adsorption  system  and  the  next  two  are  the  relaxed  energies  of 
the  isolated  molecule  and  cluster.  AEbsse  is  the  correction  for 
basis  set  superposition  error  (BSSE),  which  is  obtained  via  the 
usual  counterpoise  method.  In  general,  no  correction  is  applied 
for  zero-point  vibrational  energy  (ZPE);  however,  an  estimate 
of  the  magnitude  of  AZ?zpe  will  be  given  below.  A  negative 
AHads  constitutes  an  exothermic  process. 


3.  Construction  of  the  Model 

3.1.  The  a-SiOi  Model.  The  a-SiCF  lattice  has  been  modeled 
by  Van  Ginhoven  et  al.26  who  showed  that  the  physical 
properties  of  the  bulk  glass  can  be  accurately  described  by 
averaging  over  a  finite  set  of  three-dimensionally  periodic 
structures  formed  from  properly  constructed  small  cells.  Each 
cell  is  formed  by  cooling  (in  a  molecular  dynamics  (MD) 
calculation)  a  gas  of  24  Si  and  48  O  atoms  under  control  of  an 
empirical  potential  energy  function.  Constraints  are  imposed 
only  on  the  atom  density  and  on  the  minimum  Si— O,  O—O, 
and  Si— Si  nearest-neighbor  distances.  Because  the  process  is 
random,  each  72  atom  SF4O48  cell  is  unique  and  represents  an 
arbitrarily  chosen  small  piece  of  the  bulk  glass.  In  modeling 
bulk  glass,26  each  of  ten  different  cells  synthesized  in  this 
manner  is  then  used  to  form  a  “crystal”  with  periodic  boundary 
conditions  and  subjected  to  further  annealing  using  a  DFT  MD 
treatment.  The  average  structural  properties  (radial  distribution 
function,  ring  size  distribution,  etc.)  of  these  ten  different  glasses 
are  in  excellent  agreement  with  experiment  and  with  empirical 
MD  simulations  for  much  larger  samples  (493  Si02  units). 

The  present  work  uses  one  of  the  72  atom  clusters  as  the 
model  system  (cf.  Figures  2  and  3)  from  which  are  cut  small 
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Figure  3.  Similar  to  Figure  2  but  showing  a  reactive  region  with  a 
relatively  low  OH  density.  Diagrams  (a)  and  (b)  show  large  and  small 
models  for  the  reactive  region,  respectively.  The  asterisk  in  (a)  marks 
the  OH  site  used  for  TCP  adsorption  (see  text).  Si  atoms  that  are 
common  to  all  three  diagrams  are  shown  in  green.  In  (c),  the  view  is 
from  a  perspective  different  from  that  in  Figure  2b. 

sections  to  represent  the  reactive  adsorption  sites.  A  similar 
approach  has  already  been  used27  to  treat  the  self-trapped  exciton 
in  SiCK  The  model  is  essentially  a  cluster  of  finite  size  that,  to 
the  extent  possible,  exhibits  structural  characteristics  expected 
for  an  a-Si02  surface.  This  includes  a  distribution  of  different 
Si— OH  configurations  as  well  as  siloxane  rings  of  varying  size 
and  shape.  Hence,  the  model  permits  a  comparison  of  the 
adsorption  characteristics  of  several  different  types  of  surface 
sites  with  no  alteration  of  the  model  itself.  The  present  approach 
differs  from  that  of  other  studies  in  which  a-Si02  is  modeled 
as  a  small  cluster  containing  only  a  few  Si04  tetrahedra 
(sometimes  only  one)  or  as  a  crystalline  structure  with  a  single, 
well-defined  siloxane  ring  geometry.  A  similar  structure,28  an 
OH-free  Si24048  cluster,  has  been  used  to  treat  the  reaction  of 
H2O  with  defects  on  the  a-SiCB  surface.  Here,  the  model  was 
formed  by  cutting  an  Si24048  cluster  from  crystalline  a-quartz 
and  then  subjecting  it  to  simulated  annealing. 

In  the  present  case,  all  dangling  Si  (O)  bonds  on  the  periphery 
of  the  cluster  were  saturated  with  OH  (H).  Dangling 
— O— Si(OH)3  units,  which  are  known29  to  be  unstable  on  the 
a-Si02  surface  above  ~300  °C,  were  removed  and  replaced  with 
—OH.  The  resulting  cluster  then  has  a  stoichiometry  of 
Si2i056H28-  Next,  the  positions  of  the  O  and  H  atoms  in  the 
OH  groups  were  optimized  at  either  the  semiempirical  PM3 
level30  or  the  ab  initio  restricted  Hartree— Fock  (RHF)/3-21G 
level  while  the  rest  of  the  cluster  remained  fixed.  Two  reactive 
sites  were  then  selected.  One  (Figure  2a)  is  a  site  of  relatively 
high  OH  density  and  consists  of  an  isolated  (— O— )3Si— OH 
with  adjacent  geminal  (— O— )2Si(OHb  groups.  The  other 
(Figure  3a, b)  is  a  site  of  relatively  low  OH  density  and  consists 
mainly  of  isolated  Si— OH  groups.  Many  other  combinations 
of  isolated  and  geminal  groups  are  possible,  but  these  two  have 
been  chosen  as  representative  of  the  hydroxylated  ^-SiCL  surface 


(see  below).  A  fully  dehydroxylated  n-Si02  surface,  with  only 
(—  O— )3Si — O— Si(— O— )3  bridges,  does  not  occur  naturally  and 
requires  vacuum  annealing3132  at  ~1000  °C  to  form. 

The  OH  density  on  a  fully  hydroxylated  a-SiOi  surface  has 
been  found  experimentally32  to  be  ~5  OH  nm-2.  Heating  in 
vacuo  to  ~450  °C  is  a  typical  substrate  pretreatment  in  IR 
studies  of  TCP  and  DMMP  adsorption  (see  below).  This  desorbs 
any  physisorbed  H2O  and  converts  most  H-bonded  Si— OH 
groups  to  Si— O— Si  bridges  via  H2O  elimination,  leaving  a 
coverage  of  non-H-bonded  groups32  equal  to  about  2  OH  nm-2. 
For  the  model  used  here,  the  total  coverage  is  about  4.0  OH 
nm-2,  based  on  the  surface  area  of  a  sphere  enclosing  the  cluster. 
The  cluster  consists  predominantly  of  isolated  and  geminal 
Si— OH  groups,  and  only  two  pairs  of  Si— OH  groups  are 
sufficiently  close  to  H  bond  to  each  other.  The  0— H  bond  length 
between  two  H-bonded  silanol  sites  is  taken  here  to  be  2.0  A, 
based  on  results18  at  the  MP2/6-31G(d)  level  for  H20  and 
CH3OH  interacting  with  Si(OH)4.  Such  sites  are  not  used  in 
the  present  work  because  the  available  experimental  data  for 
TCP  and  DMMP  adsorption  (see  below)  all  focus  on  a-SiOi 
surfaces  with  a  low  coverage  of  H-bonded  OH  groups.  If  the 
H-bonded  sites  are  excluded,  the  coverage  on  the  model  cluster 
becomes  about  3.4  OH  nm-2.  Hence,  the  model  corresponds 
fairly  closely  to  the  surfaces  used  in  experimental  studies, 
although  the  coverage  of  non-H-bonded  Si— OH  groups  is  ~50% 
higher  than  on  the  real  material. 

3.2.  ONIOM  Treatment.  In  the  next  step,  an  ONIOM 
procedure33  is  applied  in  which  the  small  reactive-site  subcluster 
(Figure  2a  or  3a, b)  and  the  adsorbate  define  the  “model  system” 
while  the  entire  cluster  plus  adsorbate  comprises  the  “real 
system”.  During  geometry  optimization,  all  atoms  in  the 
ONIOM  model  system  are  allowed  to  relax  without  constraint 
while  the  rest  of  the  cluster  remains  fixed  in  position  to  prevent 
distortion  of  the  cluster  from  its  bulk-glass  configuration.  The 
ONIOM  energy  is  then  given  by  E  =  E'™  +  £^el  -  £^"del. 
“Low”  refers  to  a  low-level  treatment  which,  in  the  present  work, 
is  either  semiempirical  (PM330  or  MNDO34)  or  ab  initio  (RHF/ 
3-21G).  “High”  refers  to  a  high-level  ab  initio  treatment  which, 
in  the  present  work,  is  mainly  at  the  B3LYP/6-31+G(2d,2p)  level 
with  an  “ultrafine”  integration  grid  and  “tight”  convergence 
criteria  for  geometry  optimization  to  increase  the  level  of 
accuracy  in  the  treatment  of  weak  H-bonding  interactions.  Thus, 
in  the  definition  of  AZ?ads  given  above,  Lfmol  +  cluster)  is  the 
ONIOM  energy  with  the  adsorbate  present,  and  ii(cluster)  is 
the  ONIOM  energy  without  the  adsorbate.  The  effects  of  the 
choice  of  high-  and  low-level  method  will  be  discussed  below. 
Because  of  program  limitations,  AZTbsse  cannot  be  obtained 
directly  for  an  ONIOM  model.  Hence,  after  optimizing  the 
ONIOM  adsorption  structure  the  model  system  is  detached  and 
H  atoms  used  to  saturate  the  dangling  Si— O  bonds.  The 
counterpoise  procedure  is  then  applied  to  this  isolated  structure 
to  get  AIibsse-  Hydrogen  link  atoms  were  used  in  calculations 
for  the  ONIOM  model  system. 

Solans-Monfort  et  al.35  have  compared  ONIOM  and  periodic- 
lattice  calculations  for  the  adsorption  of  NH3  or  H2O  at  acidic 
OH  sites  in  the  aluminosilicate  material  “chabazite”.  They 
concluded  that  the  two  approaches  give  comparable  geometries 
but  that  use  of  RHF/3-21G,  rather  than  a  semiempirical  method, 
as  the  ONIOM  low-level  method  is  needed  to  obtain  A/:il(js  in 
agreement  with  the  periodic— lattice  result.  Furthermore,  it  was 
found  that  good  ONIOM  results  could  be  obtained  by  optimizing 
the  geometry  using  MNDO  as  the  low-level  method  followed 
by  a  single-point  calculation  with  RHF/3-21G  as  the  low-level 
method.  Such  single-point  calculations  of  A Z?ads  converged 
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rapidly  to  the  periodic— lattice  result  as  the  size  of  the  ONIOM 
model  system  increased.  As  a  low-level  method,  MNDO 
appeared  to  be  somewhat  superior  to  AMI.  It  should  be  noted 
that  a  periodic— lattice  approach  might  not  be  practicable  for 
the  adsorbates  of  interest  here.  Most  CWA  molecules  are  even 
larger  than  Sarin  (Figure  1),  and  very  large  unit  cells  would  be 
required  to  prevent  interaction  between  nearest-neighbor  ad¬ 
sorbates. 

Conclusions  similar  to  those  mentioned  above  were  reached 
by  Roggero  et  al.20  who  compared  ONIOM  results  for  a 
“cagelike”  cluster  model  of  Si02  with  those  obtained  in  a  fully 
ab  initio,  high-level  treatment  of  the  entire  cluster.  In  particular, 
it  was  found  that  ONIOM  results  with  RHF/3-21G  as  the  low- 
level  method  were  in  excellent  agreement  with  those  for  the 
fully  high-level  treatment  and  that  the  ONIOM  results  converged 
rapidly  as  the  size  of  the  model  system  increased.  It  was  also 
found20  that  PM3  performs  relatively  poorly  as  an  ONIOM  low- 
level  method  for  Si02-  However,  Michalkova  et  al.3  found  that 
PM3  and  RHF/3-21G  give  similar  results  when  used  as  the 
ONIOM  low-level  method  in  modeling  adsorption  of  Sarin  on 
the  aluminosilicate  material  “dickite”.  The  difference  may  be 
related  to  the  size  of  the  ONIOM  model  system,  which  was 
quite  large  in  ref  3  versus  here  and  in  ref  20.  It  is  likely  that 
with  a  sufficiently  large  model  system  the  accuracy  of  the  low- 
level  method  becomes  relatively  unimportant.  A  comparison  of 
AMI,  PM3,  and  MNDO/d  in  predicting  the  structures  of  small 
Si-containing  molecules  has  also  been  given  by  Bartlett  et  al.36 

In  the  present  work,  a  calculation  for  the  bare  (i.e.,  adsorbate- 
free)  cluster  (Figure  2a)  with  PM3  as  the  ONIOM  low-level 
method  and  B3LYP/6-31+G(2d,2p)  for  the  high-level  gave  a 
Mulliken  charge  of  g™°del  =  +1.28  |e|  for  Si  atoms  in  the 
model  system.  The  same  calculation  with  MNDO  (RHF/3-21G) 
as  the  low-level  method  gave  g™°del  =  +1.83  |e|  (+2.03  |e|). 
The  RHF/3-21G  result  is  virtually  identical  to  the  Mulliken 
charge  of  +2.036  |e|  found  in  a  periodic  RHF  treatment37  of 
bulk  a-Si02-  Hence,  c/™kIlI  is  sensitive  to  the  choice  of  low- 
level  method,  and  even  a  simple  ab  initio  method  gives  a  charge 
result  close  to  that  obtained  for  a  bulk  lattice.  Further  discussion 
is  given  below  regarding  the  different  results  obtained  for  the 
three  low-level  methods. 

4.  Results  and  Discussion 

4.1.  TCP  Adsorbed  on  a-SiCh.  The  main  interest  in  the 
present  work  is  in  the  adsorption  of  phosphonates  on  the 
hydroxylated  ci-SKT  surface.  To  our  knowledge  there  are  few 
previous  theoretical  studies14-19  of  this  particular  problem.  The 
adsorption  of  TCP  was  taken  as  a  starting  point  because  it  is  a 
smaller  and  simpler  molecule  than  either  DMMP  or  Sarin; 
hence,  test  calculations  can  be  done  more  quickly.  Furthermore, 
the  mode  of  TCP  adsorption  (P=0— H— O  bonding)  is  unam¬ 
biguous.38-39  Thus,  it  can  be  used  to  probe  differences  in  the 
reactivity  of  various  regions  of  the  model  cluster  more  easily 
than  can  a  larger,  polyfunctional  molecule  like  DMMP. 

The  goal  of  the  TCP  studies  is  to  assess  the  computational 
procedures  before  proceeding  to  the  larger  and  more  complex 
adsorbates.  Experimental  data,  in  the  form  of  IR  vibrational 
spectra,  have  been  reported38-39  for  TCP  adsorbed  on  high- 
surface-area  (HSA)  «-Si02  powders.  The  substrate  preparation40 
used  in  these  studies  (vacuum  annealing  at  400—450  °C)  yields 
a  material  with  predominantly  non-H-bonded  Si— OH  groups 
that  may  be  either  isolated  or  geminal.  Structural  data,  obtained 
from  microwave  spectroscopy,41  are  also  available  for  gas-phase 
TCP.  These  give  r(P=0)  =  1.455  ±  0.005  A  and  r(P-Cl)  = 
1.989  ±  0.002  A  for  the  bond  lengths  and  a  Z(C1— P— Cl)  bond 
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Figure  4.  Schematic  diagram  showing  the  structure  of  Cl3P=0 
adsorbed  on  the  Sbt+Hg  cluster.  The  configuration  shown  results  from 
geometry  optimization  at  the  B3LYP/6-31+G(2d.2p)  level.  The  heavy 
green  line  shows  the  chemisorption  bond. 

angle  of  103.7  ±  0.2°.  In  the  present  work,  the  corresponding 
values  computed  at  the  B3LYP/6-31+G(2d,2p)  level  are  1.460 
and  2.036  A  and  103.6°. 

4.1.1.  Small-Cluster  Model,  As  a  further  simplification,  the 
present  study  began  with  the  SisOvHg  cagelike  cluster 
(Figure  4)  described  by  Civalleri  et  al.16  Among  the  different 
such  clusters  considered  by  this  group  this  was  the  smallest  that 
gave  results  for  adsorption  of  NH3,  which  were  reasonably 
independent  of  cluster  size.  Table  1  shows  that  there  is  little 
difference  among  the  A£^ds  values  for  different  basis  sets. 
Both  the  TCP  and  the  entire  cluster  were  allowed  to  relax  during 
optimization.  From  the  computed  vibrational  modes  (at  the 
6-3 1  +G(d,  p)  level)  of  the  adsorbate  +  cluster,  bare  cluster,  and 
free  adsorbate,  AZ?zpe  =  +0.81  kcal/mol  was  obtained  after 
scaling  the  mode  frequencies  by  a  factor  of  0.9632  (see  above), 
resulting  in  A B^ds  =  —3.7  kcal/mol.  A£zpe  is  not  expected  to 
depend  significantly  on  the  basis  set.  Thus,  the  corresponding 
result  for  6-31+G(2d,2p)  was  A£Zpe  =  +0.70  kcal/mol,  giving 
AH^ds  =  — 3.6  kcal/mol.  Hence,  the  A£^ds  results  remain 
nearly  identical  after  further  correction  for  AZ?Zpe-  Figure  4 
shows  the  optimized  structure  for  adsorbed  TCP. 

A  calculation  was  also  done  at  the  6-31+G(2d,2p)  level  for 
adsorption  by  P— Cl- H— O  bonding.  The  result  was  AZ?ads  = 
—0.46  kcal/mol  and  r(C\— H)  =  3.490  A,  indicating  only  a  weak 
interaction.  This  is  consistent  with  the  smaller  charge  on  Cl 
versus  O,  as  shown  by  the  atomic  polarizability  tensor42  (APT) 
results  for  the  free  molecule  in  Table  2.  A/:a(|s  for  this  type  of 
bonding  corresponds  to  ~10%  of  the  adsorption  energy  due  to 
P=0— H— O  interaction.  The  possibility  of  adsorption  via  an 
electrostatic  interaction  between  the  highly  charged  P  atom  (cf. 
Table  2)  and  the  O  of  the  Si— OH  group  was  also  investigated, 
but  an  even  smaller  interaction  was  found  (AZ?ads  ^  —0.23  kcal/ 
mol,  r(P— O)  ~  4.389  A).  Furthermore,  a  structure  of  this  type 
would  be  sterically  hindered  for  the  larger  DMMP  and  Sarin 
molecules  (cf.  Figure  1). 

Table  1  shows  that  proper  treatment  of  anharmonicity 
(specifically,  of  the  change  in  VeXe  caused  by  adsorption)  is  a 
significant  factor  in  Av(0— H),  as  noted  in  previous  stud¬ 
ies. 15-16-21-22  The  results  for  the  bare  surface,  v0i(O— H)  and 
do  not  depend  strongly  on  the  choice  of  basis  sets,  but  An(P= 
O)  and  Anoi(0— H)  are  in  somewhat  better  agreement  with 
experiment  for  6-31+G(2d,2p).  The  voi(0— H)  and  ryce  results 
are  close  to  those  reported16  by  Civalleri  et  al.  (3762  and  78 
cm-1,  respectively)  for  the  SisOvHs  cluster  at  the  B3LYP/ 
Dunning  DZP(6d)  level  and  with  experimental  results38-39-43  for 
isolated  OH  groups  on  HSA  «-Si(+  powders.  After  adsorption 
of  TCP,  a  VeXe  of  about  106  cm-1  is  computed  for  the  three 
basis  sets. 

In  comparing  observed  and  calculated  results  for  Anoi(0— 
H),  it  must  be  recalled  that  the  former  are  obtained  at  room 
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TABLE  1:  DFT/B3LYP  Results  for  TCP  Adsorption  on  the  SisCLHs  Cluster  for  Different  Basis  Sets" 

method  AEaisb  r(P=0 — H)  Ar(P=0)c  Avas(PCl3)c  v0i(O-H)''  tye'  Av0i(O-H/ 


6-31+G(d,  p) 

-5.9  (-4.5) 

1.867 

-15 

+  13 

3740  (3759) 

76 

-223  (-168) 

6-311++G(d,  p) 

-6.4  (-4.1) 

1.880 

-16 

+8 

3746  (3786) 

74 

-221  (-179) 

6-31+G(2d,2p) 

-4.9  (-4.3) 

1.897 

-20 

+  13 

3751  (3769) 

77 

-239  (-181) 

Experiment* 

-36 

0 

3747 

-251 

Experiment'’ 

-33 

+  11 

3747 

-247 

Experiment' 

3745 

85 

“  AEads  is  in  kcal/mol,  r(P=0 — H)  is  the  H-bond  distance  in  A,  and  v,  v^xe,  and  Ar  are  in  cnT1.  ''Values  for  A£^ds  (the  BSSE-corrected  A,Eal]S) 
are  in  parentheses.  cAv(P=0)  and  Avas(PCl3)  are  the  frequency  shifts  (relative  to  the  gas  phase)  obtained  for  purely  harmonic  potentials  after 
scaling  as  described  in  the  text.  <Voi(0— H)  is  the  anharmonic  SiO— H  frequency  for  the  bare  cluster  as  described  in  the  text.  The  value  in  parentheses 
is  the  frequency  obtained  for  a  purely  harmonic  potential  after  scaling.  '’Anharmonicity  parameter  for  the  bare  cluster.  After  adsorption,  v^xe  is  about 
106  cm-1  for  the  different  basis  sets.  'At’oiiO— H)  is  the  change  in  Voi(0— H)  caused  by  adsorption  of  TCP.  Numbers  in  parentheses  are  the 
corresponding  results  obtained  after  scaling  for  a  purely  harmonic  potential.  *Data  for  bare  a-SiO:  and  for  adsorbed  TCP  from  ref  38  (recorded  at 
room  temperature).  ''Data  for  bare  a-SiO:  and  for  adsorbed  TCP  from  ref  39  (recorded  at  RT).  'Data  for  bare  a-SiCF  from  ref  43  (recorded  at  RT). 
VeXe  =  (2voi-Vo2)/2  was  obtained  from  the  reported  voi(0— H)  and  Vo2(0— H)  frequencies. 


TABLE  2:  APT  Atomic  Charges  in  Free  Molecules" 

molecule  P*=0  P=0*  R— O*  P— X* 

C13P=0  +2.10  -0.76  -0.45 

DMMP  +2.12  -0.82  -1.00 

Sarin  +2.13  -0.84  -1.09  -0.67 

"All  values  are  obtained  at  the  B3LYP/6-31+G(2d,2p)  level  with 
an  ultrafine  integration  grid.  The  asterisk  marks  the  atom  being 
described.  R  =  CHj-  for  DMMP  and  (CH,)2CH-  for  Sarin.  X  =  Cl 
for  TCP  and  F  for  Sarin. 

temperature  while  the  later  pertain  to  T  =  0.  It  has  been  shown 
experimentally44  for  NH3  adsorbed  on  HSA  a-SiCL  powders 
that  Ai'(0— H)  =  —700  cm-1  at  room  temperature  (RT)  and 
—900  cm-1  at  4  °K.  With  large  basis  sets,  a  periodic— lattice 
model  and  the  inclusion  of  anharmonicity,  the  DFT/B3LYP 
result15  for  AvoKO— H)  is  —733  cm-1.  After  various  corrections 
are  applied,  this  becomes  —633  cm-1,  which  is  still  close  to 
the  experimental  RT  value  but  underestimates  the  T  =  0  value. 
No  IR  data  are  available  for  TCP  adsorbed  on  SiCL  at  cryogenic 
temperatures,  but  the  results  in  ref  44  suggest  that  the  apparent 
good  agreement  between  theory  and  experiment  for 
Avoi(0— H)  in  Table  1  may  be  misleading.  Nevertheless,  the 
agreement  is  better  when  anharmonicity  is  explicitly  included. 

To  determine  whether  explicit  treatment  of  anharmonicity 
would  also  improve  the  agreement  between  the  observed  and 
calculated  Av(P=0),  a  similar  analysis  was  done  for  this  mode. 
Examination  of  the  normal-mode  displacements  for  TCP,  either 
free  or  adsorbed  on  the  SisCpHs  cluster,  shows  that  the  P=0 
stretch  is  essentially  decoupled  from  other  modes;  hence,  the 
approach  outlined  above  is  applicable.  In  this  case,  the  anhar¬ 
monicity  is  much  smaller  than  for  the  O— H  stretch.  For  free 
TCP,  i'eXe  =  7.0  cm-1  for  v(P=0),  using  6-31+G(2d,2p)  basis 
sets,  vs  77  cm-1  (Table  1)  for  v(0— H)  of  the  bare  Si— OH 
group.  The  approximation,  noted  above,  of  describing  the 
i'(P=0)  normal  mode  entirely  in  terms  of  displacement  of  the 
lighter  atom  is  less  valid  here  than  in  the  case  of  v(0~ H). 
Nevertheless,  the  computed  anharmonic  Voi(P=0)  for  free  TCP, 
1268  cm-1,  is  closer  to  the  experimental  gas-phase  value38,39 
of  1322  cm-1  than  is  the  result  (v(P=0)  =  1246  cm-1)  of 
scaling  the  purely  harmonic  frequency.  However,  the  adsorption- 
induced  shift  obtained  with  explicit  inclusion  of  anharmonicity, 
At'oi(P=0)  =  — 17  cm-1,  is  negligibly  different  from  the  value 
of  —20  cm-1  (Table  1)  obtained  simply  by  scaling  the  purely 
harmonic  frequency  shift.  Thus  the  effect  of  adsorption  on  vexe 
appears  to  be  much  smaller  for  v(P=0)  than  for  v(0— H). 

Finally,  Table  1  shows  that  all  three  models  give  reasonable 
results  for  Avas(PCl3),  the  shift  relative  to  the  gas  phase  for  the 
asymmetric  PCh  stretching  mode.  This  mode  is  doubly  degen¬ 
erate  in  the  gas  phase  but  splits  by  3  to  6  cm- 1  in  the  adsorbed 


state  (depending  on  the  method  of  calculation)  when  the  C3„ 
symmetry  of  the  free  molecule  is  eliminated.  The  shifts  given 
in  Table  1  are  based  on  the  average  of  the  two  frequencies  for 
the  adsorbed  species. 

4.1.2.  Large-Cluster  Model.  The  small  cluster  described  above 
may  appear  adequate  for  describing  the  adsorption  of  TCP; 
however,  deficiencies  in  this  model  will  be  identified  below. 
In  any  case  larger,  polyfunctional  molecules  such  as  DMMP 
and  Sarin  require  a  model  with  multiple  OH  groups.  Even  for 
a  relatively  simple  molecule  like  TCP,  a  realistic  treatment  of 
adsorption  must  allow  the  possibility  of  interaction  with  adjacent 
sites  not  directly  bonded  to  the  adsorbate.35  Several  tests  of  the 
adsorption  of  TCP  on  the  large-cluster  model  will  now  be 
performed.  These  include  a  comparison  of  the  high-OH-density 
(Figure  2a)  and  low-OH-density  (Figure  3a,b)  models,  a 
comparison  of  isolated  and  geminal  sites,  and  an  assessment  of 
the  effects  of  the  size  of  the  ONIOM  model  system. 

Table  3  shows  results  for  TCP  adsorption  at  the  isolated 
Si— OH  group  of  the  model  shown  in  Figure  2a.  These  were 
obtained  after  optimization  using  semiempirical  low-level 
methods.  Consistent  results  for  AE^s  and  r(P=0— H— O)  were 
found  only  with  a  subsequent  single-point  calculation  with  RHF/ 
3-21G  for  the  low-level  method.  Without  the  single-point 
calculation,  MNDO  gives  a  shorter  bond  length  but  a  smaller 
AEads  than  does  PM3.  In  both  cases,  the  optimized  structure 
(not  shown)  gave  no  indication  of  interaction  with  sites  other 
than  the  adsorption  site.  For  example,  the  nearest  adjacent  OH 
site  is  located  such  that  the  resulting  OH— 0=P  distance  is  3.98 
A,  and  any  other  relevant  distances  (e.g.,  OH— Cl)  are  even 
greater. 

In  view  of  the  results  in  Table  3,  the  9™°del  values  noted 
above  and  the  discussions  in  refs  20  and  35,  all  subsequent 
ONIOM  work  will  make  use  of  RHF/3-21G  as  the  low-level 
method.  One  adjustment  was  necessary,  namely,  the  use  of  a 
3-21G(d)  basis  set  for  the  P  atom.  In  the  molecules  of  interest 
here,  P  is  hypervalent,  having  one  ji-  and  three  cr-bonds.  This 
requires  a  higher  degree  of  variational  freedom  (i.e.,  d-orbitals) 
in  the  P  basis  set,  without  which  SCF  convergence  problems 
were  often  encountered.  In  the  following  (and  in  Table  3), 
3-21G(d)  for  the  P  atom  should  be  understood  whenever 
reference  is  made  to  3-2 1G  basis  sets  for  the  low-level  method. 

Table  4  and  Figure  5  show  results  for  TCP  adsorption  at 
different  sites,  which  indicate  some  significant  points.  First  is 
the  fact  that  for  adsorption  at  the  isolated  Si— OH  site  of 
Figure  2a,  A Eads  for  optimization  with  RHF/3-21G  as  the 
low-level  method  (—9.5  kcal/mol)  is  close  to  the  single-point 
result  after  optimization  with  MNDO  as  the  low-level  method 
(—10.1  kcal/mol,  Table  3).  This  is  consistent  with  previous 
results35  for  a  different  system  (NH3  adsorption  on  an  alumi- 
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TABLE  3:  Results  for  TCP  Adsorption  on  Large-Cluster  a-SiCL  Models  for  Different  ONIOM  Methods'1 


method 


r(P=0--H— O) 


B3LYP/6-3 1  +G(2d,2p)  :PM3 
B3LYP/6-31+G(2d,2p):PM3//:RHF/3-21G‘' 
B3LYP/6-31+G(2d,2p):MNDO 
B3LYP/6-31+G(2d,2p):MNDO//:RHF/3-21G 


-10.9  (-10.5)  2.001 

-2.9  (-2.5) 

-5.8  (-5.3)  1.918 

-10.1  (-9.7) 


“  Energies  are  in  kcal/mol;  bond  lengths  are  in  A.  In  all  cases,  adsorption  occurs  at  the  isolated  OH  of  the  model  shown  in  Figure  2a.  '’Numbers 
in  parentheses  are  the  BSSE-corrected  (A£^ds)  values.  The  same  A£bsse  (+0.42  kcal/mol)  is  used  in  all  cases  because  the  high-level  method  is  the 
same.  The  notation  means  that  the  geometry  was  first  optimized  using  B3LYP/6-3 1+G(2d,2p)  for  the  high-level  and  PM3  for  the  low-level  method. 
Then  a  single-point  calculation  was  done  with  RHF/3-21G  as  the  low-level  method  (with  3-21G(d)  for  the  P  atom). 


TABLE  4:  Results  for  TCP  Adsorption  on  Different  a-SiCL 
Large-Cluster  ONIOM  Models" 

model4  A  £adsc  r(P=0— H-O) 


(i)  isol  OH  (Figure  2a)  -9.5  (-9.0) 

(ii)  gem  OH  (Figure  2a)  —6.2  (—5.7) 

(iii)  isol  OH;  large  model  (Figure  3a)  —  1 1 .0  (—10.4) 

(iv)  isol  OH;  small  model  (Figure  3b)  —7.6  (—7.0) 


1.845 

1.877 

2.148,  2.203d 
1.804 


“  All  calculations  are  done  at  the  B3LYP/6-31+G(2d,2p):RFIF/3-21G 
level  (with  a  3-21G(d)  basis  set  for  the  P  atom).  Energies  are  in  kcal/ 
mol,  and  bond  lengths  are  in  A.  4isol  =  isolated;  gem  =  geminal.  The 
structures  are  labeled  (i),  etc.,  to  facilitate  referencing  in  the  text  and 
in  Figure  5.  ^Values  in  parentheses  are  A£^ds  for  AEbsse  = 
+0.42  kcal/mol.  This  structure  exhibits  H  bonding  between  P=0  and 
two  Si— OH  groups. 


Figure  5.  Optimized  structure  for  TCP  adsorbed  at  the  Si— OH  sites 
indicated  in  Table  4.  The  labels  (i),  etc.,  correspond  to  those  in  Table 
4.  For  clarity,  only  the  ONIOM  model  systems  are  shown,  and  each  is 
in  approximately  the  same  orientation  as  in  Figure  2a  or  3a,b.  The 
additional  H  atoms  used  as  ONIOM  "link”  atoms  are  not  shown.  The 
heavy  green  lines  show  the  dominant  bonding  interactions. 

nosilicate  material)  and  serves  as  a  useful  test  of  the  present 
computational  procedure. 

In  all  cases  shown  in  Table  4,  A£^ds  is  larger  than  the  value 
of  —4.3  kcal/mol  found  for  the  small  SisCFHg  cluster 
(Table  1).  This  may  be  because  of  the  rigidity  of  the  SisCLHs 
structure,  which  consists  of  fused  siloxane  rings  with  only  three 
or  four  Si  atoms,  as  opposed  to  the  “loose”  and  open  structure 
of  the  large-cluster  model.  This  permits  only  limited  relaxation 
in  response  to  adsorption,  even  though  the  entire  SisCLHs  cluster 
is  unconstrained  during  optimization.  In  this  context,  one  notes 
that  of  structures  (i),  (ii),  and  (iv)  the  most  stable  is  structure 


(i)  in  which  adsorption  occurs  at  the  center  of  the  fully  relaxed 
subcluster.  When  only  one  Si— OH  group  is  involved  in 
H-bonding,  geminal-  and  isolated-OH  sites  have  been  found45 
to  give  similar  AE^s  values.  Hence,  the  difference  in  AEads 
between  (i)  and  (ii)  is  not  thought  to  be  due  to  the  presence  of 
the  second,  noninteracting  geminal  OH  group  in  (ii). 

The  optimized  structures  in  Figure  5  exhibit  differing  degrees 
of  interaction  between  TCP  and  multiple  OH  sites.  This 
emphasizes  a  basic  difficulty  in  treating  a  disordered  system 
with  a  small  finite-cluster  model.  For  example,  (iii)  shows  a 
“bidentate”  structure  with  H  bonding  between  TCP  and  two 
Si— OH  sites.  This  can  be  achieved  only  for  a  local  geometry 
in  which  two  such  sites  are  present  with  the  correct  separation 
and  relative  orientation.  Note  that  adsorption  at  the  geminal  site 

(ii)  does  not  lead  to  this  energetically  most  favorable  bidentate 
structure.  Other  structures  suggest  a  possible  contribution  from 
O— H— Cl— P  interaction  (not  shown)  involving  OH  sites  outside 
of  the  ONIOM  model  system,  although  this  should  be  relatively 
weak,  as  noted  above.  One  structure,  (ii),  shows  no  indication 
of  such  interactions.  None  of  these  can  be  taken  to  represent  a 
unique  and  universal  model  for  TCP  adsorption  on  ci-SiCF. 

Taken  together,  the  results  in  Tables  1  and  4  suggest  that 
there  is  a  requirement  for  a  sufficiently  flexible  structure,  one 
with  adequate  freedom  to  relax  in  the  vicinity  of  the  adsorption 
site  and  with  an  ONIOM  model  system  that  includes  all  sites 
that  might  be  involved  in  adsorbate  bonding.  When  these 
conditions  are  reasonably  well-satisfied,  there  is  not  a  great 
difference  in  AEads  between  structures  with  quite  different 
geometries  (e.g.,  (i)  and  (iii)  in  Figure  5).  A  more  precise 
treatment  of  TCP  adsorption  would  focus  on  a  range  of  different 
adsorption  sites,  all  modeled  in  accord  with  the  criteria  stated 
above.  This  would  presumably  decrease  the  spread  in  AEa(js 
values  shown  in  Table  4.  However,  the  main  interest  here  is  in 
the  adsorption  of  DMMP  and  Sarin.  Further  refinements  in  the 
treatment  of  TCP  adsorption  are  therefore  not  justified  but  the 
results  obtained  above  will  be  applied  to  modeling  the  adsorbates 
of  interest. 

Finally,  the  SiO— H  stretching  frequency  was  examined  for 
the  case  of  adsorption  with  the  structure  in  Figure  5(i).  For  the 
bare  cluster,  voi(0— H)  =  3744  and  VeXe  =  77  cm-1  were 
obtained  which  are  in  good  agreement  with  experiment 
(Table  1).  Adsorption  of  TCP  gave  Avoi(0— H)  =  —267  cm-1. 
A  similar  calculation  for  structure  (ii),  which  has  the  smallest 
AEads  of  those  shown  in  Figure  5,  gave  Avoi(0— H)  = 
—260  cm-1.  These  values  (which  pertain  to  T  =  0)  are  slightly 
larger  than  the  RT  experimental  result  of  about  —250  cm-1 
(Table  1).  Bearing  in  mind  the  remarks  made  previously 
regarding  the  temperature  dependence  of  Afoi(0— H),  the 
agreement  with  experiment  is  reasonable. 

It  is  useful  to  compare  the  results  for  AEads,  r(P=0 — H— O), 
and  Anoi(0— H)  in  Table  1  with  those  for  structures  (i),  (ii), 
and  (iv)  given  above  and  in  Table  4.  In  all  cases, 
/-(P=0— H— O)  and  Avoi(0— H),  which  directly  reflect  the 
strength  of  the  H-bond  interaction,  are  similar,  whereas,  the 
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AZiads  values  differ.  For  molecules  in  solution,  the  Badger— 
Bauer  relationship12  predicts  a  strong  and  nearly  linear  increase 
in  Avoi(0— H)  with  H-bond  strength.  The  present  results  suggest 
that  AZiads  depends  on  other  factors  in  addition  to  the  H-bond 
strength,  such  as  relaxation  at  the  adsorption  site  and  possibly 
interaction  of  the  adsorbate  with  adjacent  sites.  Thus,  the  small 
SisCbHg  cluster  performs  fairly  well  in  describing  the  strength 
of  a  single  H  bond  but  not  in  predicting  the  adsorption  energy. 

4.2.  DMMP  Adsorbed  on  a-SiCh.  The  insights  gained  in 
the  study  of  TCP  adsorption  will  now  be  applied  to  the  case  of 
DMMP.  The  approach  will  be  to  use  the  large-cluster  model  to 
determine  A£ads,  the  adsorption  geometry  and  Aroi(0— H)  and 
to  use  the  SisCbHg  cluster  in  computing  the  internal  vibrational 
modes  of  adsorbed  DMMP.  The  latter  step  is  taken  for  two 
reasons.  First,  even  with  the  use  of  analytical  gradients,  a 
complete  normal-mode  calculation  for  the  large-cluster  model 
is  too  computationally  expensive.  Second,  the  partial  geometry 
optimization  performed  for  the  large-cluster  model  (see  above) 
precludes  an  accurate  calculation  of  the  lower-energy  modes 
of  adsorbed  DMMP  because,  unlike  the  SiO— H  stretch,  many 
of  these  are  coupled  to  displacements  of  atoms  in  the  cluster. 
The  foregoing  discussion  of  TCP  indicates  that  reasonable 
results  are  to  be  expected  for  the  internal  modes  of  an  adsorbate 
bonded  to  the  SisOvHs  cluster.  Experimental  and  theoretical 
results  for  the  conformation  and  vibrational  modes  of  gas-phase 
DMMP  have  been  discussed  in  ref  10.  As  in  ref  10,  the  gas- 
phase  DMMP  conformer  (Figure  lb)  used  in  the  present  work 
is  the  one  identified  previously46  as  the  lowest  in  energy. 

4.2.1.  Experimental  Results.  The  interaction  of  DMMP  with 
a-SiCB  HSA  powder  has  been  studied  by  Kanan  and  Tripp38 
using  IR  spectroscopy  and  by  Henderson  et  al.47  using  temper¬ 
ature-programmed  desorption  and  Auger  electron  spectroscopy. 
The  IR  results  are  interpreted  in  terms  of  non-dissociative 
adsorption  via  H-bond  formation  between  isolated  or  geminal 
Si— OH  sites  and  the  two  molecular  CH3O—  groups.  The  P=0 
group  is  thought  not  to  be  directly  involved  in  adsorption.  The 
activation  energy  for  desorption  is  estimated47  to  be  ~  16.9  kcal / 
mol. 

The  main  effects  in  the  IR  spectrum  that  result  from 
adsorption  are  the  following: 

(a)  a  shift  of  the  SiO— H  stretch  by  Av(0— H)  =  —524  cm-1, 
relative  to  the  bare  surface; 

(b)  a  shift  of  the  SiO— H  bend  by  A<3(0— H)  ~  +160  cm-1, 
relative  to  the  bare  surface; 

(c)  a  shift  of  the  P=0  stretch  by  Av(  P=() )  =  — 19  cm-1, 
relative  to  the  gas-phase; 

(d)  little  or  no  change  in  the  va(C— H)  and  vs(C— H) 
antisymmetric  and  symmetric  stretches  of  the  —  OCH3  groups; 

(e)  a  weakening  and/or  broadening  of  the  vs  and  14  symmetric 
and  antisymmetric  P— O— CH3  stretching  modes,  relative  to  the 
gas  phase,  to  the  point  of  disappearance. 

As  a  further  point,  it  is  noted  that  H  bonding  between  DMMP 
and  H2O  is  found48  both  theoretically  and  experimentally  to 
occur  via  a  P=0— H— O  linkage  with  experimental  values  of 
Av(P=0)  =  — 17  cm-1  and  Av(0— H)  =  —203  cm-1  for  the 
symmetric  stretching  mode  of  H20.  Studies14-19  at  the  B3LYP/ 
6-31+G(d,  p)  level  of  the  interaction  between  orthosilicic  acid, 
Si(OH)4,  and  the  dimethyl  phosphonate  ion,  [(CHsO^PC+l-, 
also  indicate  an  energetic  preference  for  H  bonding  to  the  two 
lone  O  atoms  versus  to  the  methoxy  O  atoms.  In  IR  spectro¬ 
scopic  studies  of  the  interaction  of  DMMP  with  organic  self- 
assembled  monolayers  terminated  in  OH  tail  groups49  and  with 
polymeric  siloxane  materials,50  bonding  is  suggested  to  occur 


Figure  6.  Initial  and  final  (optimized)  structures  for  adsorbed  DMMP. 
For  clarity,  only  the  molecule  and  the  ONIOM  model  system  (cf. 
Figure  3a)  are  shown.  The  additional  H  atoms  used  as  ONIOM  “link” 
atoms  are  not  shown.  Diagrams  (a)  and  (c)  show  the  initial  structures, 
and  (b)  and  (d)  show  the  optimized  structures.  The  heavy  green  lines 
show  the  dominant  bonding  interactions. 

through  P=0— H— O  formation,  which  leads  to  shifts  in  v(P= 
O)  as  large  as  —48  cm-1  relative  to  the  gas  phase. 

Thus  the  interpretation  proposed  by  Kanan  and  Tripp,38 
involving  adsorption  exclusively  via  the  DMMP  methoxy 
groups,  differs  from  models  put  forth  in  connection  with  other 
similar  systems.  However,  strong  support  for  this  interpretation 
comes  from  IR  data38  for  the  adsorption  of  (CHsYSilOCHsY-* 
( x  =  1,  2,  3).  The  modes  involving  — OCH3  groups  behave 
identically  to  those  in  DMMP,  although  Av(0— H)  is  only 
—397  cm-1  versus  —524  cm-1  for  DMMP.  In  the  case  of 
DMMP,  only  one  shifted  SiO— H  stretch  is  observed,  not  two 
as  would  be  expected  if  adsorption  involved  both  a  strong  bond 
to  CH3O—  and  a  weak  bond  to  P=0.  It  therefore  appears 
that  the  surface  chemistry  of  a-Si02  might  differ  from  that  of 
the  model  systems  discussed  above. 

4.2.2.  Computed  Results.  Calculations  were  done  starting  with 
a  model  in  which  adsorption  occurs  via  a  single  P=0— H— O 
bond  (Figure  6a)  or  a  single  (CH3)(P)0— H— O  bond 
(Figure  6c).  As  before,  B3LYP/6-31+G(2d,2p)  was  used  for  the 
ONIOM  high-level  method,  and  RHF/3-21G  [3-21G(d)  for  P] 
was  used  for  the  low-level  method.  The  a-Si02  model  in  all 
cases  is  that  shown  in  Figure  3a.  The  results  are  summarized 
in  Table  5.  In  the  first  case,  rearrangement  of  the  initial  structure 
occurs  to  give  a  final  optimized  structure  (Figure  6b)  with 
bonding  of  the  P=0  to  two  Si— OH  groups,  as  did  TCP 
(Figure  5(iii)).  The  results  obtained  are  AZ?ads  =  —21.6  kcal/ 
mol  (AZ4ds  =  —20.0  kcal/mol)  and  1.916  and  1.936  A  for  the 
two  H-bond  lengths.  In  the  second  case,  the  optimized  structure 
(Figure  6d)  shows  a  short  H  bond  (1.783  A)  between  one 
CH3O—  and  an  Si— OH  group  and  a  longer  bond  (2.134  A) 
between  the  other  CH3O—  and  a  second  Si— OH.  This  is 
essentially  the  structure  proposed  by  Kanan  and  Tripp.38  The 
adsorption  energy  in  this  case  is  AZ?ads  =  —15.1  kcal/mol  (A 
Zs^ds  =  —13.6  kcal/mol).  These  AZs^ds  values  are  both  larger 
than  the  corresponding  result  for  TCP  (—10.4  kcal/mol. 
Table  4)  which  is  consistent  with  the  experimental  observation38 
that  DMMP  adsorbs  more  strongly  on  n-Si02  than  does  TCP. 
The  A/vds  values  are  also  close  to  the  observed  desorption 


Adsorption  of  TCP,  DMMP,  and  Sarin  on  Amorphous  Si02 

TABLE  5:  Results  for  DMMP  Adsorption  on  a-SiOj" 

AFC  b 
^^ads 

r(0— H)*-' 

Avoi(O-H)' 

■>  A 6(0— FI)'  Av(P=0)' 

Figure  6b  —20.0 

1.916 

1.936 

-420 

~+257  -29 

Figure  6d  —13.6 

1.783 

2.134 

-423 

-+170  +7 

experiment/ 

-524 

~+160  -19 

“  Energies  are  in  kcal/mol,  bond  lengths  are  in  A,  and  frequencies 
are  in  cm-1.  ''Values  obtained  for  adsorption  on  the  large-cluster  model 
shown  in  Figure  3a.  'Two  H  bonds  are  formed  in  either  case  (see  text). 
The  two  H-bond  lengths  for  each  structure  are  tabulated.  ''Values 
obtained  from  an  exact  treatment  of  anharmonicity  (see  text)  for 
adsorption  on  the  large-cluster  model  shown  in  Figure  3a.  Only  the 
stronger  (i.e.,  shorter)  of  the  two  H  bonds  is  considered.  'Values 
obtained  for  the  small-cluster  ShCFHg  model.  The  calculated  harmonic 
values  have  been  scaled  by  a  factor  of  0.9632  (see  text).  'Data  from 
ref  38. 

energy47  of  about  16.9  kcal/mol,  and  including  a  small  A£zpe 
term  (<  1  kcal/mol,  see  above)  would  improve  the  agreement 
in  the  case  of  adsorption  via  the  P=()  group. 

A  calculation  was  also  done  starting  with  a  structure  (not 
shown)  in  which  DMMP  adsorbs  via  short  H  bonds  (~1.8  A) 
between  both  CH3O—  groups  and  two  isolated  Si— OH  groups. 
This  was  done  in  an  effort  to  see  if  a  more  stable  form  of  the 
structure  shown  in  Figure  6d  could  be  forced  by  starting  closer 
to  the  desired  final  configuration.  Early  in  the  geometry 
optimization,  the  DMMP  rearranged  from  the  initial  structure 
to  one  exhibiting  an  H  bond  to  a  single  CH30—  group  and  a 
second  bond  between  another  Si— OH  and  the  P— O  group.  The 
resulting  A£ads  was  —10.1  kcal/mol,  indicating  a  metastable  (i.e., 
local  minimum)  configuration  relative  to  those  shown  in 
Figure  6b,d.  In  effect,  the  H  bond  to  either  group  (P=0  or 
CH3O— )  interferes  with  optimum  bonding  to  the  other.  It  thus 
appears  that  the  structure  shown  in  Figure  6b  is  the  most  stable 
that  can  be  achieved. 

The  relevant  vibrational  properties  can  now  be  studied 
(Table  5).  Anharmonic  Avoi(0— H)  values  for  the  structures 
shown  in  Figure  6b, d  were  obtained  as  described  above.  In  either 
case,  only  the  shorter  (stronger)  H  bond  was  considered,  and 
this  gave  Ai’oi(0— H)  =  —420  cm-1  (—423  cm-1)  for  the 
structure  in  Figure  6b  (Figure  6d).  These  are  to  be  compared 
with  the  experimental  value  of  Ai'oi(C)— H)  =  —524  cm-1.  As 
is  typically  found,1516’21  the  calculated  values  (which  are  nearly 
identical  for  the  two  structures)  underestimate  somewhat  the 
experimental  result. 

The  internal  vibrational  modes  of  DMMP  were  computed  for 
adsorption  in  either  of  the  two  possible  configurations  at  the 
isolated  Si— OH  group  of  the  SisOyHs  cluster  (as  discussed 
above).  The  calculation  was  done  at  the  B3LYP/ 
6-31+G(2d,2p)  level  with  complete  relaxation  of  both  cluster 
and  adsorbate.  The  AZ?ads  values  were  —9.7  and  —4.8  kcal/mol, 
respectively,  for  H  bonding  to  P=0  and  to  CH3O— .  The  H-bond 
lengths  were  1.700  and  1.853  A,  respectively.  As  was  found 
for  TCP  adsorption,  the  small-cluster  energies  are  significantly 
smaller  in  magnitude  than  those  for  the  larger  and  more  flexible 
cluster  with  two  H  bonds  per  DMMP.  Nevertheless,  DMMP 
adsorbed  via  the  P=0  group  is  again  found  to  be  the  more 
favorable  of  the  two  structures,  which  is  consistent  with  the 
large-cluster  results  given  above.  A  i'(  P=() )  for  H  bonding  to 
the  P=0  group  is  —29  cm-1,  whereas  for  H  bonding  to  the 
CH3O—  group  it  is  +7  cm-1.  The  experimental  result  is 
Ae(P=0)  =  — 19  cm-1.  Thus,  H  bonding  to  the  CH3O—  group 
gives  a  Ae(P=0)  of  the  wrong  sign. 

For  gas-phase  DMMP,  the  out-of-phase  and  in-phase  v(C~ 
O)  modes  are  computed  (at  the  B3LYP/6-31+G(2d,2p)  level) 
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to  lie  at  1007  and  1031  cm-1,  respectively,  after  scaling  as 
described  above.  The  corresponding  experimental  values49  are 
1050  and  1075  cm-1.  The  absolute  values  differ  somewhat  from 
experiment,  due  to  systematic  errors24  in  the  frequency  calcula¬ 
tion,  but  the  difference  between  the  in-phase  and  out-of-phase 
modes  is  practically  identical  in  theory  (24  cm-1)  and  experi¬ 
ment  (25  cm-1). 

For  DMMP  adsorbed  via  the  P=0  group,  these  modes  show 
little  or  no  shift  in  frequency  (~3  cm-1).  For  adsorption  via 
the  CH3O—  group,  v(C~ O)  is  at  980  cm-1  for  the  H-bonded 
CH3O—  and  1025  cm-1  for  the  free  CH30— .  This  shift  of 
45  cm-1  is  not  by  itself  expected  to  make  the  H-bonded  C— O 
stretch  undetectable  in  the  IR  spectrum.  For  either  type  of 
adsorption,  examination  of  the  normal  mode  atomic  displace¬ 
ments  shows  that  the  C— O  stretching  modes  are  coupled  to 
Si— O— Si  bond-stretching  “lattice  modes”  of  the  cluster  and/or 
to  the  Si— O— H  bending  mode,  all  of  which  occur  in  about  the 
same  energy  region.  The  coupling  can  be  described  as  “strong” 
in  the  sense  that  the  C— O  stretching  normal  modes  also  involve 
large  Si— O  stretching  and/or  Si— O— H  bending  motions.  This 
may  account  for  the  broadening  seen  experimentally38  for  the 
C— O  stretching  modes  when  either  DMMP  or  ( C  H+S  i  ( Of '  H  i  )^-x 
is  adsorbed  on  a-SiOi.  Thus  the  broadening  of  the  C— O 
stretching  modes  is  not  necessarily  an  indication  of  H  bonding 
to  the  O  atom  of  the  CH3O—  groups  in  DMMP. 

Table  5  also  gives  results  for  the  effect  of  adsorption  on 
<5(0— H),  the  frequency  of  the  Si— O— H  bending  mode.  An 
exact  value  for  A(5(0— H)  is  difficult  to  determine  in  the 
computational  results  because  of  strong  coupling  between  this 
mode  and  the  Si— O  and  C— O  stretching  modes,  as  discussed 
above.  For  the  bare  cluster,  (5(0— H)  is  found  at  832  cm-1.  For 
DMMP  adsorbed  via  the  P=0  group.  Si— O— H  bending 
contributes  significantly  to  several  modes  in  the  range  of  about 
1002  to  1089  cm-1.  The  highest  of  these  appears  to  be  almost 
pure  Si— O— H  bending.  Using  this  value  results  in  an  estimate 
of  A(5(0— H)  ~  +257  cm-1.  In  a  sense,  the  mixing  of 
Si— O— H  bending  with  the  C— O  and  Si— O  stretching  modes 
is  the  result  of  adsorption,  which  shifts  the  bending  mode  into 
the  energy  range  of  the  stretches.  Similar  considerations  for  the 
case  of  adsorption  via  the  CH3O—  group  give  A(5(0— H)  ~ 
+  170  cm-1,  based  on  the  frequency  of  the  mode  that  most 
closely  approximates  a  pure  Si— O— H  bending  motion.  Although 
this  result  is  very  close  to  the  experimental  value  (~+160  cm-1), 
the  results  for  A£^ds  and  for  Av(P=0)  continue  to  argue  in 
favor  of  H  bonding  to  the  P=0  group. 

4.3.  Sarin  Adsorbed  on  n-SiOi.  There  are,  to  our  knowledge, 
no  experimental  studies  of  Sarin  adsorbed  on  a-SKK  Michalk- 
ova  et  al.9  have  reported  theoretical  results  for  the  interaction 
of  Sarin  with  Si(OH)4.  Bonding  occurs  via  O— H— 0=P 
interaction  (with  an  H-bond  distance  of  1.756  A)  accompanied 
by  an  additional  interaction  between  an  H  of  the  —  CH3  group 
and  the  O  of  a  second  Si— OH  group  (at  a  distance  of  2.657  A). 
H  bonding  to  the  O  atom  of  the  iPr—  O— P  group  is  not  found 
to  be  important.  At  the  B3LYP/6-31G(d)  level,  A£^ds  = 
—7.6  kcal/mol  is  obtained.  Experimental  and  theoretical  results 
for  the  conformation  and  vibrational  modes  of  free  Sarin  have 
been  discussed  in  ref  10.  The  conformer  used  for  gas-phase 
Sarin,  shown  in  Figure  lc,  is  identical  to  that  used  in  ref  10. 
This  is  the  second-lowest-energy  conformer,51’52  which  lies  <0.2 
Kcal/mol  above  the  absolute  lowest-energy  conformer.  Different 
Sarin  conformers  are  defined  in  terms  of  rotations  about  the 
iPr— O—  and  — O— P  bonds,  and  different  enantiomers  of  the 
same  conformer,  obtained  by  interchanging  the  F  and  phos- 
phonyl  O  atoms,  are  isoenergetic.51 
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TABLE  6:  Results  for  Sarin  Adsorption  on  n-SiOj" 

H-bond  site  A £ads  A£^ds  r(H— X)4 

P=0  -21.4  -19.9  1.779,2.133 

iPr-O-P-F  -11.3  -9.5  1.944  (O),  2.146  (F) 

P-F  -12.4  -11.2  2.198,2.413 

a  Energies  are  in  kcal/mol,  and  bond  lengths  arein  A.  The  atom  shown 
in  bold  is  the  one  involved  in  H  bonding.  4X=0  or  F.  For  P=0  and 
P— F,  H  bonds  are  formed  to  two  Si— OH  groups.  For  iPr— O— P— F, 
one  H  bond  forms  to  the  O  atom  and  another  to  the  F  atom. 


Figure  7.  Similar  to  Figure  6b  but  showing  the  optimized  structure 
for  Sarin  adsorbed  via  H  bonding  to  the  P=0  group.  The  violet  sphere 
is  the  F  atom. 

In  the  present  work,  three  modes  of  adsorption  were 
considered.  These  are  H  bonding  to  the  O  atom  of  the  P=0  or 
iPr— O— P  group,  as  for  DMMP,  and  also  to  the  F  atom.  In  each 
case,  the  large-cluster  model  (Figure  3a)  is  used,  and  the  starting 
structure  involves  only  one  H  bond  (as  in  Figure  6a,c).  The 
results  for  the  final  optimized  structures  are  summarized  in 
Table  6.  As  in  the  case  of  TCP  and  DMMP,  the  most  stable 
structure  (Figure  7)  involves  H  bonding  between  the  P=0  group 
and  two  Si— OH  groups.  Adsorption  via  the  F  atom  also  involves 
H  bonding  between  the  F  and  two  Si— OH  groups  but  is  less 
stable  than  bonding  to  P=0.  The  optimized  structures  were 
examined  for  the  possibility  of  H  bonds  to  other  functional 
groups;  however,  all  such  interatomic  distances  were  3.2  A  or 
greater.  In  comparison  to  the  H-bond  distances  shown  in 
Tables  1,  4,  and  6,  this  is  considered  to  be  too  large  for  a 
significant  contribution  to  AZ?ads.  The  initial  structure  with  a 
single  H  bond  to  the  O  atom  of  the  iPr— O— P  group  rearranged 
during  optimization  to  one  with  this  H  bond  intact  and  with  a 
second  H  bond  between  F  and  another  Si— OH.  However,  this 
is  again  less  stable  than  two  H  bonds  to  the  P=0  group.  The 
vibrational  properties,  computed  for  the  most  stable  optimized 
structure  using  the  procedure  described  above,  were 
Avoi(0— H)  =  —510  cm-1  (for  the  shorter  of  the  two  H  bonds) 
and  Ai'(P=0)  =  —25  cm-1.  Given  the  results  for  DMMP  (Table 
5),  the  computed  Avoi(0— H)  for  Sarin  probably  underestimates 
the  experimental  value,  and  the  computed  Ai'(P=0)  probably 
overestimates  the  experimental  value. 

5.  Summary 

Ab  initio  QC  calculations  have  been  performed  to  study  the 
interaction  of  TCP,  DMMP,  and  Sarin  via  H  bonding  to  Si— 
OH  groups  on  the  hydroxylated  a-SiCL  surface.  The  results  are 
as  follows. 

(1)  A  small,  cagelike  cluster  gives  reasonable  results  for  the 
vibrational  spectrum  of  adsorbed  TCP  and  DMMP.  This 
includes  the  adsorption-induced  red  shift  of  the  silanol 
v(0— H)  stretching  mode  when  anharmonicity  is  treated  explic¬ 


itly.  However,  the  adsorption  energy  is  underestimated  because 
of  the  rigidity  of  the  model,  which  limits  the  extent  of  relaxation 
at  the  adsorption  site,  and  to  the  omission  of  interactions  between 
the  adsorbate  and  more  than  one  Si— OH. 

(2)  For  the  adsorbates  studied  here,  an  ONIOM  treatment 
requires  the  use  of  an  ab  initio  method  (here  RHF/3-21G)  for 
the  low  level,  rather  than  a  semiempirical  method.  Including  d 
orbitals  in  the  low-level  basis  set  (e.g.,  3-21G(d))  is  required 
for  the  adsorbate  P  atom. 

(3)  For  all  three  species  studied  here,  H-bond  formation 
between  the  O  atom  of  the  P=0  group  and  two  Si— OH  groups 
is  the  most  energetically  favorable  adsorption  mechanism.  For 
DMMP  and  Sarin,  H  bonding  to  alkoxy  O  atoms  is  less 
favorable,  as  is  H  bonding  to  the  Sarin  F  atom. 

(4)  The  energetics  and  geometry  of  adsorption  for  DMMP 
and  Sarin  on  a-SiO?  depend  on  the  local  Si— OH  environment. 
The  bidentate  P=0(— H— 0)2  structure  is  the  dominant  factor 
in  adsorption.  However,  contributions  from  other,  weaker 
interactions  can  be  expected  in  real  systems  if  the  local 
Si— OH  configuration  is  such  that  these  effects  are  energetically 
feasible.  H  bonding  to  more  than  one  functional  group  is  seen 
in  the  present  work  for  adsorption  of  DMMP  in  a  metastable 
configuration  with  one  H  bond  to  the  O  atom  of  a  CH3— O— P 
group  and  a  second  to  the  P=0.  It  is  also  seen  for  Sarin 
adsorption  (again  in  a  metastable  configuration)  via  the  O  atom 
of  the  iPr— O— P  group,  where  a  second  H  bond  to  the  F  atom 
also  forms.  Similar  effects  have  also  been  observed  in  compu¬ 
tational  studies3,4,9  of  Sarin  on  other  oxides. 

(5)  The  characteristics  of  DMMP  and  Sarin  with  respect  to 
adsorption  on  the  a-SiCF  surface  are  closely  similar.  In 
particular,  the  computed  adsorption  energies  of  the  two  species 
are  virtually  identical.  Hence,  DMMP  is  a  good  simulant  for 
Sarin  in  this  regard.  It  is  not,  however,  to  be  concluded  that  the 
same  is  necessarily  true  for  all  adsorbent  materials  and  under 
all  conditions.  For  example,  the  effect  of  coadsorbed  H2O  on 
the  adsorption  of  DMMP  or  Sarin  has  not  yet  been  analyzed, 
either  computationally  or  experimentally.  Furthermore,  for 
adsorption  on  OH-free  y-A^Ch,  which  occurs  via  a  different 
mechanism,10  DMMP  may  interact  somewhat  more  strongly 
than  does  Sarin. 
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